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Utilizing an eigenfunction decomposition, we study the growth and spectra of energy in the vortical (geostrophic) and wave 
(ageostrophic) modes of a three-dimensional (3D) rotating stratified fluid as a function of e = f/N, where / is the Coriolis param- 
eter and N is the Brunt- Vaisala frequency. Throughout we employ a random large scale forcing in a unit aspect ratio domain and set 
these parameters such that the Froude and Rossby numbers are roughly comparable and much less than unity. Working in regimes 
characterized by moderate Burger numbers, i.e. Bu = 1/e 2 < 1 or Bu > 1, our results indicate profound change in the character of 
vortical and wave mode interactions with respect to Bu = 1. Previous analytical work concerning the qualitatively different nature of 
these interactions has been in limiting conditions of rotation or stratification domination (i.e. when Bu < 1 or Bu 2> 1 respectively). As 
with the reference state of e = 1, for e < 1 the wave mode energy saturates quite quickly and the ensuing forward cascade continues to 
act as an efficient means of dissipating ageostrophic energy. Further, these saturated spectra steepen as e decreases: we see a shift from 
A: -1 to fc -5 / 3 scaling for kf < k < k a (where kf and k^ are the forcing and dissipation scales, respectively). On the other hand, when 
e > 1 the wave mode energy never saturates and comes to dominate the total energy in the system. In fact, in a sense the wave modes 
behave in an asymmetric manner about e = 1. With regard to the vortical modes, for e < 1, the signatures of 3D quasigcostrophy are 
clearly evident. Specifically, we see a A: -3 scaling for kf < k < k^ and, in accord with an inverse transfer of energy, the vortical mode 
energy never saturates but rather increases for all k < kf. In contrast, for e > 1 and increasing, the vortical modes contain a progressively 
smaller fraction of the total energy indicating that the 3D quasigcostrophic subsystem plays an energetically smaller role in the overall 
dynamics. Combining the vortical and wave modes, the total energy for k > kf and e < 1 shows a transition as k increases wherein the 
vortical modes contain a large portion of the energy at large scales, while the wave modes dominate at smaller scales. There is no such 
transition when e > 1 and the wave modes dominate the total energy for all k > kf . 

1 Introduction 

In terms of the Froude (Fr) and Rossby (Ro) numbers, geophysical phenomena occur in an environment 
characterized by strong rotation and stratification, i.e. these numbers are O(10 -1 ) or smaller (Gill 1982). 
In the present work, we consider a forced scenario where /, N (the Coriolis parameter and Brunt- Vaisala 
frequency) are comparable, and thus define these parameters by Fr = (efkj-) 1 / 3 /N and Ro = (efkj-) 1 / 3 / f 
where £/, are the energy input rate and forcing wavenumber, respectively. The influence of rotation and 
stratification on fluid flows has been examined in detail, and recent reviews include Riley and Lelong (2000) 
for stratification, Cambon et al. (2004) for rotation and Cambon (2001) for both stratification and rotation. 

For purely rotating flows (Q = Qz), the system supports inertial waves obeying the dispersion relation 
cr(k) = ±2Qk z /k. When forced at small scales and for Ro < 0(1), the principal effect of rotation is seen 
in the two-dimensionalization of the flow: simulations indicate the spontaneous formation of columnar 
structures from initially isotropic conditions (Smith & Waleffe 1999, see Bartello, Metais & Lesieur 
(1994) for analogous behavior in the decaying problem). In effect, there is a transfer of energy to the 
slow wave modes (k z = => °~(k) = 0) and this transfer has an inverse character. Given that exactly 
resonant interactions cannot transfer energy to slow wave modes (Greenspan 1969), near-resonances play 
an important role in this process (Smith & Lee 2005). Analytically, for Ro — > 0, the limiting dynamics 
involving purely resonant interactions contain three parts (Babin, Mahalov & Nicolaenko 1999). The first 
part, consisting of zero frequency interactions, yields the two-dimensional (2D) Navier-Stokes equations. 
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The second part, involving interactions of the fast-zero-fast type, results in the passive driving of the 
vertical velocity by means of the aforementioned 2D flow (Babin, Mahalov & Nicolaenko 1996, Chen et 
al. 2005). Finally, including the third category of interactions involving three fast wave modes results 
in the so-called 2^ dimensional equations (Babin, Mahalov & Nicolaenko 1999). When forced at large 
scales (for Ro < 0(1)), in agreement with the analytical picture of so-called "phase turbulence" for 
fast wave modes with a ^ (Babin, Mahalov & Nicolaenko 1996), the transfer of energy to smaller 
scales is seen to be inefficient (Cambon et al. 1997, see also Belief et al. 2006 for statistical closures of 
these wave-turbulence phenomena). Also, for rapid rotation, along with an accumulation of energy near 
the forcing scale there is phenomenological (Zhou 1995) and numerical support for a k~ 2 scaling in an 
intermediate range of scales (Yeung & Zhou 1998, Belief et. al. 2006). 

For uniform vertical stratification, in addition to gravity waves satisfying cr(k) = ±Nkh/k (where 
k\ = k 2 x + ky) the equations possess a distinct zero frequency mode, the so-called vortical mode (Lelong & 
Riley 1991). With forcing at small scales and for Fr < 0(1), the principal effect of stratification is seen in 
the layering of fields: numerical simulations indicate the spontaneous emergence of layers perpendicular to 
the ambient stratification (Metais & Herring 1989, Smith & Waleffe 2002; see Godeferd & Cambon 1994, 
Riley & De Bruyn Kops 2003 for the 3D decaying case; see Smith 2001 and Sukhatme & Smith 2007b for, 
respectively, the forced 2D case and the decaying 2D case). As for pure rotation, there is an an inverse 
transfer of energy to the slow wave modes (k^ = => decoupled layers) (Smith & Waleffe 2002). In the 
limit Fr — » 0, the presence of a reduced system consisting of diffusively coupled 2D layers is anticipated 
by a scaling analysis of the original equations (Riley, Metcalfe and Weissman 1981, Lilly 1983, Majda & 
Grote 1997). With forcing at large scales and Fr < 0(1), there are significant differences for the scaling 
of the energy spectra depending on the nature and dimensionality of the force (Waite & Bartello 2004, 
Lindborg 2006). However, in accord with the analytical picture of "unfrozen" wave mode cascades (Babin, 
Mahalov & Nicolaenko 1997), simulations employing these vastly differing forcing mechanisms agree upon 
the robust nature of energy transfer to small scales. 

In the presence of both rotation and stratification, the focus of the present work, an appropriate 
governing framework is provided by the 3D rotating Boussinesq equations. Note that the aforementioned 
results regarding rotating flows hold when N = 0, and that it is not immediately clear how the limit of 
large e = f/N converges to a purely rotating flow. In essence the role of the extra vortical mode that 
exists whenever N ^ needs careful consideration^] On the other hand, the limit of small Fr with fixed 
Ro has been rigorously analyzed and shown to obey a splitting between the vortical and fast wave modes 
wherein both the vortical and slow wave modes evolve as vertically sheared horizontal flows (VSHF) 
(Embid & Majda 1998). Setting / = in this limiting case results in the previously mentioned small Fr 
reduced system of diffusively coupled 2D layers for purely stratified flows. 

In the linear and non-dissipative limits, the 3D rotating Boussinesq equations support wave-like motions 
that obey the dispersion relation (see Majda 2003 for a recent overview) 



where e = f /N = Fr/Ro. Here modes corresponding to o± are called wave (or ageostrophic) modes 
while the mode corresponding to do is referred to as the vortical (or geostrophic) mode (Bartello 1995, 



1 Quite interestingly, a family of dynamically active dispersive scalars (which includes the familiar barotropic beta plane equations as a 
member) is seen to have a similar structure of resonant interactions (Sukhatme & Smith 2007a). Indeed, there the analog of the present 
transfer to columnar structures is the spontaneous generation of zonal flows. 

1 From recent discussions with Prof. A.J. Majda, the limit of -Ro — > while holding Fr ~ 0(1) appears to be described by a reduced 
system involving barotropic dynamics coupled with oscillatory density and vertical velocity fields. This is presently being examined in 
greater detail. 



a±(k) = ± 



N(kj + e 2 ^) 1 / 2 
k 



a {k) = 



(1) 
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Smith & Waleffe 2002). In a unit aspect ratio domain, the Burger number is defined as Bu = (Ro/Fr) 2 . 
In the presence of strong stratification and rotation a regime relevant to atmosphere-ocean dynamics is 
Fr,Ro <C 1 with Bu ~ 1 (Gill 1982). A rigorous analysis of the limit Fr,Ro — > 0, Bu ~ 1 indicates 
a splitting between the vortical and wave modes of the 3D rotating Boussinesq system wherein the 
vortical modes evolve as an independent 3D quasigeostrophic (QG) system while the wave modes evolve 
via (catalytic) wave-vortical-wave and wave-wave-wave interactions (Babin et al. 1997, Embid & Majda 
1998). Bartello (1995) performed a detailed study of geostrophic adjustment for / ~ N with a± ~ ±iV, 
implying a system with vortical and nearly constant frequency wave modes. Examining the decay of 
a large scale initial condition, he noted the aforementioned decoupling between the vortical and wave 
modes. Further, apart from the 3D QG dominance of the vortical modes (accompanied by a k~ 3 scaling) 
it was seen numerically that catalytic wave-vortical-wave interactions dominated the wave mode dynamics 
(accompanied by a k^ 1 scaling) resulting in a rapid transfer of energy to small scales. In this case, the 
catalytic interactions are responsible for the dissipation of ageostrophic energy and the emergence of a 
geostrophically adjusted state (Bartello 1995). In a later study, Waite & Bartello (2006) found a gradual 
transition from stratified to quasigeostrophic turbulence as the Rossby number was decreased for a fixed 
small Froude number. Employing a small scale random forcing for 1/2 < e < 2, evidence for the dominance 
of vortical modes by QG-like dynamics was also presented in Smith & Waleffe (2002), who detected 
the presence of a /c~ 5 / 3 upscale energy cascade amongst the 3D vortical modes (see also Metais et al. 1996). 

As e gradually deviates from unity, the wave mode frequencies are non-zero and no longer constant, i.e. 
f < \&±\ < N (and vice versa if N < /). As is known (and will be elaborated on in the next section), 
the dependence of wave mode frequencies on wavenumber places restrictions on the existing vortical-wave 
mode interactions while allowing for additional interactions not present when / = N. As part of an 
extensive body of work, Babin et al. and Embid & Majda have examined a variety of asymptotic limits 
for both the purely rotating and the rotating stratified Boussinesq equations; this work is reviewed in 
Babin, Mahalov and Nicolaenko (2002) and Majda (2003) respectively. In particular, rigorous analysis 
with respect to the Bu <C 1 and Bu 3> 1 limits (note that the relevant Froude and/or Rossby numbers 
are always assumed to be small) indicates a profound difference in the behavior of the vortical and wave 
modes in these rotation vs. stratification dominated cases, respectively. Indeed, based on these limiting 
scenarios, our aim is to examine if the difference in the behavior of the modes begins to gradually manifest 
itself as one moves in the direction e > 1 or e < 1. Of course, as Bu = 1/e 2 , a moderate deviation from 
unity in e is much more substantial in terms of Bu. As suggested by the so-called "regime diagram" 
in the review by Babin et al. (2002), even though rotation domination is analytically justifiable only 
in the limit of small Bu, it is interesting to investigate how quickly rotation dominates as the Burger 
number decreases from unity (and vice versa for the stratification dominated regime). In fact, these 
intermediate regimes in Bu are precisely those that are of geophysical interest. Further, intermediate 
Bu regimes are not presently accessible to theoretical analysis, hence motivating a numerical investigation. 

Utilizing a random, isotropic, large scale forcing that projects equally on all the modes and keeping both 
Fr,Ro <C O(l), our objective in the present work is to test if we can detect a difference in the behavior 
of the various modes with respect to e = 1, with e = 1 as a well understood reference state (Bartello 
1995). The reasoning behind employing an isotropic random force is to avoid biasing any specific mode, 
i.e. the system is allowed to choose the dominant response. Specifically, numerical simulations are used 
to investigate how varying e (such that we go from Bu = 1/25 to Bu = 25) affects the scaling of energy 
spectra, the growth of energy and the overall role played by the vortical and wave modes of the 3D rotating 
Boussinesq system. Further, this lack of severe anisotropy enables us to confidently use a numerical grid 
with uniform spacing in all directions. Indeed, if one proceeds to extremes in e it is quite important to 
consider skewed grids that are able to properly resolve the dynamics (see for example Lindborg 2006 for the 
strongly stratified case). In the next section we introduce the basic equations and their limiting form for 
Ro, Fr — > as an eigenfunction decomposition. From the latter, one anticipates a change in the dynamics 
with e. In Section III we present a systematic numerical examination as e deviates from unity. For a fair 
comparison, all simulations are carried out with a fixed input energy rate and varying /, N so as to keep 
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both Ro,Fr <C 0(1)- Finally in Section IV we collect and discuss the results. 

2 The basic equations and their spectral decomposition 

The rotating Boussinesq equations in a 3D periodic setting read 



Du 

— + fz xu = -V<t> - N6z + uV 2 u 
D0 » 



Nw = kV z 9 (2) 

Dt y ' 

where u = (u, v, w) is the 3D velocity field, (j) = p/po is a scaled pressure and / is the Coriolis parameter with 
rotation assumed to be aligned with the z-axis. Equations © result from considering periodic perturbations 
to a state of hydrostatic balance wherein the density profile satisfies p = po — bz + p' with \bz\ <C p - 
Further, we have set p' = (Np /g)6 where N = (gb/ po) 1 ^ 2 ■ Finally, v and k are the viscosity and diffusivity 
respectively. The solution to the general linear initial value problem can be represented as (see Embid & 
Majda 1996, 1998 for a detailed spectral analysis of the linear operator arising from ([2])) 

= 5^ a/3(k,t)<l>i3(k)exp[ik'x\ (3) 

k /3=0,± 

where ^ = (u v w 8) T , ao, a± are determined by the initial/forcing conditions and <j>$, <f>± are the mutually 
orthogonal eigenfunctions corresponding to <7 , a± (the forms of the eigenfunctions can be found in Embid 
& Majda 1999 and Smith & Waleffe 2002, see also Bartello 1995 for an alternate construction). As the 
basis formed by these eigenfunctions is complete, we use it to expand the solution to the nonlinear problem 
in ([2]) (Embid & Majda 1998, Smith & Waleffe 2002). Symbolically the nonlinear coupling yields 

+ ia a (k)a a (k) = ^Y1 C ti a ^>7® ( 4 ) 
where A represents a sum over p. , q such that k = p + q. The indices a,/3,j run over 0,+, — and 

kpq 

is the interaction co-efficient. Before proceeding we introduce some notation: triplets (•, •, •) refer to a 
particular type of triad interaction, where z refers to the vortical mode and / denotes the wave mode 
(not to be confused with the vertical coordinate and the Coriolis parameter, respectively). For example, 
(z, /, /) represents the evolution of a vortical mode via the interaction of two wave modes (i.e. the first 
entry evolves via an interaction of the latter two entries in the triplet). 



2.1 Resonant and near-resonant interactions 

Setting a a = b a exp (—ia a (k)t) where a = {0, +, — } in (j4]), the so-called limiting dynamics for Fr, Ro — > 
— due to the cancellations amongst the oscillatory phase factors — are given exclusively by purely resonant 
interactions, 



(5) 
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where Ar restricts the sum to triads with zero sum of wave frequencies as well as wavevectors (the reader 
is referred to Embid & Majda 1996, 1998 for a formal derivation of ([5]) via a two timescale analysis, see 
also Anile et al. (1993) for an overview of such techniques in the context of general nonlinear dispersive 
equations). The vortical modes evolve as 



dJ ^T 1 = E c l% Mp>o(?1 + {eg* Mp1M^ [ct±(pH ^ =0] (6) 

A 

while the wave mode evolution is given by 



A 



{C+t* &+($M$} [ff+ ^ +<7±( * ) ~ <7+( * )=0] + {c+-± b _(^ 6± (0}[<T_(p)+ CT± (9)-<T + (fc)=o] ( 7 ) 

and a equation analogous to (|7|) governs b-(k). As we are interested in moderate values of e with 
nonzero wave frequencies a± 6 [/, N], the (z,z,f) and (f,z,z) interactions do not appear in ([6]) and ([7]) 
respectively as they can never be resonant. 

It is important to note that ©, ([7]) are limiting equations and do not account for all of the energy 
transfer among the different modes in the finite Ro, Fr system. In particular, for Fr, Ro small but non- 
zero, near-resonances (or quasi-resonances) play a significant role in the transfer of energy in addition to 
exact resonances. For example, consider a discrete initial value problem. For near resonances with wave 
frequencies summing to 5 <C 1, and denoting min(i?o, Fr) by 7 <C 1, then the phase factors implicit in 
(j3J) oscillate with periods ~ 5/7. These near-resonances will contribute until t ~ j/5 and some routes of 
energy transfer that are precluded under resonance can potentially be adopted via these near-resonant in- 
teractions, hence altering predictions from purely resonant considerations (an example of which is provided 
by Watson & Buchsbaum 1996). d In a periodic domain with small-scale forcing, Smith & Lee (2005) have 
demonstrated that near-resonant triads are responsible for the spontaneous generation of asymmetric large 
scale columnar structures from small-scale fluctuations in purely rotating 3D flows. In fact, by broadening 
the resonance condition in the kinetic equation, near-resonant interactions have recently been shown to 
have important consequences in statistical theories for the evolution of random wave fields (Janssen 2003, 
Annenkov & Shrira 2006). The present work, as with the foundational references, assumes an inherent 
heirarchy wherein much of the dynamics can be understood via exact and near-resonant interactions. Of 
course, in reality there are contributions from interactions that are far from resonance. Section 2.2 discusses 
exact resonances in the case e = 1 (for which there are no near-resonances), and Section 2.3 discusses both 
exact and near-resonances for e ^ 1. In Section 3, the simulation results for the full dynamical equations 
(4) indicate that exact and near-resonances are indeed helpful to explain the evolution of energies and 
spectra associated with the dynamics for the parameter regimes considered herein. 



2.2 e — 1 : Vortical and constant frequency wave modes 

It is known that resonant interactions involving two wave and one vortical mode are catalytic in the sense 
that the vortical mode acts as a mediator of energy exchange between the wave modes but does not 



x In periodic domain, the importance of near-resonances is further highlighted by the fact that it may be impossible for a triad to satisfy 
the dispersion relation when the admissible wavenumbers are restricted to Z (Kartashova 1994) — see Connaughton et al. (2001) and 
Lvov et al. (2006) for examples. 



2 Of course, there is also the possibility of successive resonant triads leading to quartet interactions (Newell 1969). 
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change itself (Bartello 1995, Embid & Majda 1998). Keeping this mind, for the special case of for e = 1 
(i.e. a± = ±N), © and simplify to 



dbo(k, t) 
dl 



db + (k,t) 
dl 



A 



kpq 



(8) 
0) 



and we obtain Bartello's equations (see Eqn 19 in Bartello 1995). Note that, because there is no possibility 
of near-resonances when e = 1, the exchange of energy between the vortical and wave modes is likely to 
be a minimum, making the case e = 1 quite possibly optimal in the sense of dynamical splitting between 
vortical and wave modes even at non-zero Fr, Ro. 

In (|9"j), the vortical modes evolve as an independent 3D QG system, while the wave modes are advected 
via this QG flow (see Smith & Waleffe 2002 for an explicit calculation showing that the {z, z, z) interactions 
lead to the 3D QG potential vorticity (PV) equation). Given that the vortical modes follow the 3D QG 
PV equation, we expect a large fraction of the energy entering these modes to flow upscale (Charney 
1971). Further, the equations governing the wave modes are the same as those for a passive scalar and 
we have a situation wherein these uncoupled wave modes are advected by a predominantly large scale 
QG flow. Further, this situation where the bulk of the energy of the advecting flow is at large scales 
qualitatively corresponds to the Batchelor regime for passive scalars (Batchelor 1959, Kraichnan 1974). 
This picture was verified by Bartello's (1995) decaying simulations. Starting from a (fairly) random large 
scale initial condition (characterized by a scale I) and setting / pa N in ([2]), it was seen that the energy 
in the vortical (geostrophic) modes moved upscale. For scales smaller than /, there emerged a transient 
£;~ 3 potential enstrophy transfer spectrum. Simultaneously, the energy in the wave (ageostrophic) modes 
moved downscale to be dissipated and resulted in a transient /c _1 (passive scalar like) spectrum. 



2.3 e y£ 1 : Vortical and dispersive wave modes 

As the wave mode frequencies begin to depend on wavenumber, not only do other classes of resonant 
interactions become permissible (see below) in the evolution of the wave and vortical modes, the 
aforementioned near-resonant interactions also enter the picture. Indeed, when Fr, Ro are small (though 
^ 0) and unequal, the vortical and wave modes are expected to evolve by and large independently, but we 
do expect some transfer of energy from the vortical to wave modes and vice versa due to the non-catalytic 
nature of the (z, f, /) and (/, z, f) near-resonances. 

Vortical modes: It is expected that as long as both Fr, Ro are small and comparable the vortical mode 
dynamics will be close to quasigeostrophy — as mentioned, these are the essential ingredients required 
for the rigorous demonstration of 3D QG dominance of the vortical modes (Babin et al. 1997, Embid 
& Majda 1998). But note that in the Fr — > 0, fixed Ro limit (i.e. stratification dominated regime) the 
vortical modes (along with the = slow wave modes) evolve as a vertically sheared horizontal flow 
(Embid & Majda 1998). In contrast for Fr = oo, small Ro (purely rotating regime), the vortical mode in 
fact no longer exists. Of course these are extreme cases, but they do prompt an inquiry into the behavior 
of the vortical modes as e deviates from unity. 



1 As an aside, this exchange of energy between the different modes connects to the issue of the non-existence of an invariant slow manifold 
for the rotating Boussincsq system — see for example Vanneste (2006), and the references therein, for a recent discussion of this subject. 
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Wave modes: As compared to the f = N case, we now have the added possibility of (/, /, /) resonant 
and near-resonant interactions. However, exact (/, /, /) resonances are difficult to achieve with discrete 
wavevectors given the algebraic complexity of simultaneously satisfying the triad and dispersion relations. 
In addition, they are only possible when e < 1/2 or e > 2 (Smith & Waleffe 2002). Further, for Bu ~ 0(1), 
the analysis of Babin et al. (2002) indicates that the set of (/, /, /) near resonances is very sparse. □ 
Hence for moderate values of e, it is unlikely that these additional resonant and near-resonant (/, /, /) 
interactions would a cause significant deviation from the behavior seen for e = 1 with constant wave mode 
frequencies. 

With regard to the (/, z, /) interactions, as e deviates from unity, these resonant and near-resonant 
interactions are restricted to wave modes that satisfy \k z /kh\ ~ \<lz/Qh\- Thus the (f,z,f) exact and 
near-resonances for e / 1 are not as numerous as catalytic (f,z,f) exact resonances for the case e = 1. 
In purely rotating flow, where the vortical mode no longer exists and z denotes a zero frequency wave 
mode, this restriction (which now becomes k z /\k\ ~ q z /\q\) confines the (f,z,f) interactions only to 
re-distribute energy among similar scales (Babin et al. 2002). The result is the so-called frozen cascade 
with an inefficient ageostrophic transfer of energy from large to small scales. Indeed, the inefficiency of 
energy transfer to small scales has been repeatedly noted in numerical and statistical closure studies of 
purely rotating turbulence (Cambon et al. 1997, Yeung & Zhou 1998 and Bellet et al. 2006). Even though 
the depleted transfer to small scales is seen for e 3> 1, this leaves us with the curious problem of how the 
change in the wave mode behavior arises as one eases away from e = 1. In contrast, for a stratification 
dominated regime the presence of the vortical mode (z) leads to the "unfreezing" of the ageostrophic 
forward transfer by the (/, z, f) resonant and near-resonant interactions (Babin et al. 2002) — in the 
special f = N case this unfreezing manifests itself in catalytic (/, z, /) resonances that, as discussed, lead 
to the rapid dissipation of ageostrophic energy (Bartello 1995). 



3 Numerical work 

The numerical scheme to solve ([2]) is a de-aliased pseudo-spectral method implemented in a unit aspect 
ratio domain at a resolution of 256 x 256 x 256. Time stepping is carried out by a third order Runge-Kutta 
method. We restrict the forcing to low wavenumbers (large scales) 3 < k < 4. Further, the forcing is white 
in time and is chosen from a Gaussian density function. The rate of input energy (denoted by eA in all 
cases is approximately the same, while the values of /, N are changed so as to keep both Ro,Fr <C 0(1) 
and achieve the desired e. Specifically, fixing es, we start with (Fr,Ro) = (0.01,0.01) for e = 1; for e < 1 
we choose (Fr, Ro) = (0.01, 0.02), (0.01, 0.03), (0.01, 0.04), (0.01,0.05); vice versa for e > 1. 

The only dissipation is small scale hyperdiffusivity given by (V 2 ) 8 (with v = k). Hyperdiffusivity allows 
for a range of large and intermediate scales unaffected by small- scale damping. As one expects growth 
in the vortical mode energy for low wavenumbers k < kf, the entire flow will not achieve a statistically 
steady state without large scale damping. However, since the the interaction of large scale damping with 
large scale wave modes is not well understood, we instead allow large scales to evolve freely. We terminate 
all simulations when the energy in wavenumbers k < kf is comparable to the energy in forced modes, after 
which time finite size effects are expected to influence the dynamics. 



2 Such estimates are required for investigating the issue of global regularity for the 3D rotating Boussinesq system, recently tackled by 
Babin et al. (2002). From a physical perspective, it is the forward transfer of energy that is likely to cause a breakdown in smoothness 
of the solutions and a control on these forward transfers is essential for any such investigation. 
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3.1 Wave modes 

Fig. dH) shows the wave-mode and total energies as functions of time. As is evident, for e < 1 the wave 
mode energy saturates. In marked contrast, for e > 1, the wave modes soon come to dominate the energy 
in the system. Further, comparing the fraction of energy in the wave and vortical modes, Fig. ([2j) indicates 
that the vortical mode dominance for e < 1 begins at an earlier time as stratification becomes stronger. 
On the other hand, for e > 1, dominance of the wave modes is more pronounced with increasing rotation. 

Fig. ([3]) shows the wave mode energy spectra for e < 1 after saturation (see the caption for details). 
Closely following Bartello's observations in the decaying / = N case (Bartello 1995), for e = 1 the scaling 
attains a k~ l form for kf < k < kd consistent with the passive Batchelor regime (Batchelor 1959, Kraichnan 
1974). Interestingly, as e becomes smaller, the wave mode spectra steepen and approach a A; -5 / 3 scaling by 
e = 1/5, i.e. for (Fr, Ro) = (0.01, 0.05). For e > 1, even though the wave mode energy never saturates, the 
spectra appear to attain a near-invariant shape (first panel of Fig. (jlj)), i.e. they keep shifting upwards as 
the energy increases. Note that, like pure rapid rotation (Yeung & Zhou 1998), for larger values of e there is 
some indication that a power law scaling may be emerging for an intermediate range of scales. Comparing 
the wave mode spectra for e < 1 and e > 1, it is noticeable that the peak at the forcing scale is much 
more pronounced when e > 1, i.e. rotation inhibits energy transfer out of the forcing scale. Irrespective of 
e, the energy at wavenumbers k < kf does not grow significantly indicating the universal forward nature 
of energy transfer amongst the wave modes. 



3.2 Vortical modes 

From Fig. ([I]) and the first panel of Fig. (j2J) it is evident that the increase in total energy at larger times 
for e < 1 is due to the growth in the vortical modes. As was noted above, the wave mode energy in these 
cases saturates. Further, the growth of energy at long times for e < 1 is at wavenumbers k < kf and is 
the result of an inverse transfer of energy. As mentioned above, the simulations are terminated when the 
energy in wavenumbers k < kf is comparable to the energy in the forcing scale to avoid finite size effects. 
With regard to the scaling of the vortical mode spectra, for e < 1 (first panel of Fig. ([5])) we observe a 
scaling of the form k~ 5 for kf < k < kd, consistent with the expectations from 3D QG dynamics (Charney 
1971). 

Moving to e > 1, it is evident from the second panel of Fig. (j2j) that the vortical modes contain a small 
fraction of the total energy in the system. In fact, the fraction of energy in the vortical modes decreases 
with increasing e, i.e. in an energetic sense they play a progressively smaller role in the dynamics. Again 
the energy in the vortical modes increases at all scales k < kf and we stop the simulations before finite 
size effects are expected to contribute. Considering the spectra shown in the second panel of Fig (J3J), the 
scaling k~ 3 for kf < k < k^ suggests 3D QG dynamics (though we notice the appearance of a shallower 
form at higher wavenumbers £?2 < k < kd as e increases). 



3.3 Total energy 

Combining the vortical and wave modes, Fig (jfij) shows the total energy spectra. For e < 1, we observe an 
unequal partition of energy among the wave and vortical modes at varying scales. Specifically, the vortical 
modes have a large fraction of the total energy at large scales, while the wave modes account for most 
of the energy at small scales. Given that the vortical mode spectrum scales more steeply than the wave 
modes, this results in a gradual steep-shallow transition in the total energy spectrum, as seen in the first 
two panels of Fig @. On the other hand, for e > 1, the partition is always in the favor of the wave modes 
and the total energy spectra follow the wave spectra for all k > kf. 
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4 Results and discussion 

By a series of numerical simulations we have examined the effect of varying e = f /N on the vortical and 
wave modes of the 3D rotating Boussinesq equations. Specifically, we have focused on flows characterized 
by strong rotation and stratification, i.e. both / and N are large so as keep Ro, Fr <C 1 and roughly 
comparable. In addition to the presence of the zero-frequency vortical mode, the system supports wave 
modes with frequencies cr±(k) € [f,N]. 

For e = f/N < 1 (i.e. Bu > 1, stratification stronger than rotation), the wave mode energy saturates 
quite quickly. Much like the well understood constant frequency (e = 1) case (Bartello 1995), the resulting 
forward cascade continues to act as efficient way of removing ageostrophic energy from the system. 
Simultaneously, the vortical modes exhibit a pronounced transfer of energy to large (k < kf) scales, while 
the high wavenumber energy spectra associated with the vortical modes scale as k~ 3 for kf < k < kj. In 
essence, the picture painted in Bartello (1995) involving the 3D QG dominance of the vortical modes and 
rapid adjustment via a wave mode cascade is valid for e < 1. This picture also agrees with the analytical 
work by Babin et al. (1997) and Embid & Majda (1998) regarding the splitting of of the vortical and 
wave modes for Fr ~ Ro — > wherein the vortical modes follow 3D QG dynamics while the wave mode 
cascades are "unfrozen" and result in an efficient transfer of energy to small scales. Given that the wave 
modes saturate, it is also worth noting that our results clearly indicate the energetic dominance of the 
3D QG dynamics when both Ro, Fr are small (though non-zero) and the Burger number is equal to or 
moderately larger than unity. Interestingly, we observe that restrictions on the resonant and near-resonant 
(/, z, f) interactions when e ^ 1 appear to manifest themselves in a steepening of the high-fc wave 
mode spectra from a scaling (a nonlocally dominated Batchelor regime where the fast modes are 
predominantly driven in a passive manner) to A; -5 / 3 scaling (a local 3D turbulence like behavior where 
the fast modes play a dynamically active role). Further, we notice a bias in the partition of total energy 
among the different modes at varying scales. Specifically, the vortical modes contain a large portion of the 
total energy at scales ~ kf while the wave modes account for most of the energy at smaller scales. This 
naturally introduces a gradual steep-shallow transition in the total energy spectrum. In a broader context, 
the focus on forward transfers of energy in the present work serves as a counterpart to Smith & Waleffe 
(2002) wherein small scale random forcing was employed to probe the inverse transfer of energy for 
the regimes 1/2 < e < 2 and £«1, and to elucidate the resulting 3D QG or VSHF dominance, respectively. 

Proceeding to e > 1 (i.e. Bu < 1, rotation stronger than stratification), we immediately detect 
an asymmetry in the behavior of the wave modes. Indeed, now the wave modes never saturate and 
soon dominate the entire energy in the system. In spite of this, their spectra do appear to achieve 
near-invariance, i.e. they appear to retain their shape while shifting upwards. Further, for large e there 
are signs of an emergent (steeper than fc~ 5 / 3 ) power law for an intermediate range of scales akin to the 
weak forward cascade in purely rotating turbulence (Cambon et al. 1997, Yeung & Zhou 1998). This 
situation is similar to the observations for the shallow water equations whereby switching on rotation was 
seen to inhibit forward transfer (Yuan & Hamilton 1994; see also Farge & Sadourny 1989 for remarks on 
the difficulty in achieving geostrophic adjustment in a similar scenario). As for e < 1, the vortical modes 
continue to follow 3D QG dynamics though now their role is quite small since they contain only a small 
fraction of the total energy in the system. In fact, observing that the fraction of energy in the vortical 
modes decreases with increasing e provides some hope that, at least in an energetic sense, the e>l limit 
of the 3D rotating Boussinesq system will transition smoothly to a purely rotating flow. The scaling of 
the vortical- mode spectra for e > 1 shows a /c~ 3 form for kf < k <k^ and is followed, for larger values of 
e, by a shallower form for ki < k < kj. Finally, in contrast to e < 1, now the wave modes dominate the 
total energy in the system for all k > kf. 

With regard to atmospheric phenomena, recent very high resolution studies of the rotating Boussinesq 
equations with large scale forcing and Fr <C Ro in a skewed aspect ratio domain show a similar k~ 3 scaling 
for vortical-mode spectra, k~ 5 ^ 3 scaling for wave-mode spectra, and vortical-mode energy dominance at 
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large scales (Kitamura &; Matsuda 2006). Indeed, the spectral transitions in the total energy (also in the 
potential and kinetic energies) that occur quite naturally in the Boussinesq system (see also the discussion 
in Bartello 1995) when e < 1 are reminiscent of the classic synoptic-mesoscale Nastrom-Gage spectrum 
(Nastrom &; Gage 1985). However, it should be kept in mind that the /c -5 / 3 portion of the Boussinesq 
spectra arises from wave modes whereas some observational evidence points to a vortical mode dominance 
even at small scales (Cho, Newell k, Barrick 1999). Furthermore, as clearly put forth in the recent work of 
Tulloch & Smith (2006), it is difficult to imagine a consistent theory of the midlatitude troposphere that 
does not explicitly address the potentially complicated evolution of buoyancy (or potential temperature) 
on domain boundaries (see for example Held et al. 1995 and Sukhatme k, Pierrehumbert 2002). 
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Figure 1. Wave-mode and total energies as functions of time. Upper Panel: e < 1. The dash-dot lines with symbols are the wave-mode 
energies (lower bunch) while the solid lines with symbols (upper bunch) are the total energies. The saturation of wave-mode energy in 
all cases is quite evident. Second panel: e > 1 (same notation). As is evident, the wave modes account for most of the energy in the 

system. 
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Figure 2. Ratios of wave-mode energy to total energy (solid lines with symbols) and vortical-mode energy to total energy (dash-dot 
line with symbols). For e < 1 (upper panel), initially the wave modes have a larger fraction of the total energy but after a short time, as 
the wave modes saturate, the vortical modes proceed to contain a larger fraction of the total energy in the system. In contrast for e > 1 
(lower panel), the wave modes dominate the energy budget and in fact, the fraction of energy in the vortical modes decreases with 

increasing e, i.e. as the rotation becomes stronger. 
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Figure 3. Wave mode spectra with ln(kf) ss 3.2 for e < 1. Here, Fr = 0.01, Ro = 0.01 — > 0.05. The spectra are plotted (vertically 
shifted for clarity) after the wave mode energy saturates, and are seen to change from a A: -1 to a fc -5 / 3 scaling for kj < k < kd as e 

decreases from unity. At the start of each simulation the spectra are very steep (concentrated near the forcing scale) and become 
progressively shallower. As the energy in the wave modes saturates, the spectra cease to evolve and attain the invariant shape plotted 
here. The simulations are halted when energy starts to accumulate in the large-scale vortical modes so as to avoid finite size effects. 
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Figure 4. Wave mode spectra with ln(fcj) 3.2 for e > 1. Here ijo = 0.01, Fr = 0.01 — > 0.05. The upper panel shows the time 
evolution of the spectra for e = 5. Note that even though the energy keeps growing, the spectra appear to achieve an invariant shape. 
The lower panel shows the other e > 1 cases — for clarity they have been vertically shifted. Akin to purely rotating flows, for large e 
there is some sign of a steeper than 5/3 power law emerging at an intermediate range of scales. 
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Figure 5. Vortical mode spectra. The first panel shows the spectra for all e < 1: in accord with QG dynamics, the scaling follows a 
k~ 3 form for kj < k < fcj. The second panel shows the scaling for e > 1: the spectra again follow a k~ 3 form but appear to shallow out 
for larger wavenumbers as e increases. Note that now the vortical modes contain a much smaller amount of energy than the wave 

modes (see Fig. 
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Figure 6. Total, wave and vortical mode spectra for e = 1 (upper), e = 1/4 (middle) and e = 4 (lower) (other e < 1 and e > 1 cases 
show respectively similar behavior). For e < 1, the contribution of the vortical modes to the total energy is significant at scales ~ kt, 
but decreases as we move to progressively smaller scales — indeed, this is particularly pronounced when e = 1. On the other hand for 

e = 4 (and all other e > 1 cases) the wave modes dominate V k > kf. 



